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' The response of a correlated nuclear system to an external field is discussed. 

The Bethe-Salpeter equation for the dressed vertex is solved. The kernel of the 
integral equation for the vertex is chosen consistently with the approximation for 
the self-energy. This guarantees the fulfillment of the f-sum rule for the response 
^ ' function. 
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Processes occurring in dense nuclear matter are modified by the effects of the medium. 
This happens for neutrino rates in neutron stars [H El El IH El El E! and for particle 
or photon emission in hot nuclear matter [HI El El El E] ■ In particular the effect of 
the off-shell propagation of nucleons in the medium is important for the subthreshold 
^ ■ particle production in heavy ion collisions. The role of correlations for the neutrino 
emission could be especially important for processes in hot stars. Soft particle emission 
^ \ is influenced by the multiple scattering in the medium, which is equivalent to including 
\^ \ vertex correction modifying the coupling of the external current to the dressed nucleons 
|13j . In-medium modifications of the vertex cannot be described only as modification of 
the coupling constant. The three-point function describing the modified coupling of the 
external current to the fermions in medium has a complicated analytical structure and 
depends on the incoming momenta and energies [Hj 

For the case of weak coupling of the external current or of the produced particles 
to the nucleons, the problem is equivalent to the calculation of the response function in 
the correlated medium. For the Hartree-Fock approximation the response function can 
be calculated from the random-phase approximation, which however is usually taken at 
the ring diagrams level only. The use of advanced approximations for the description 
of the many-particle system with in-medium propagators dressed by the interaction and 
scattering, raises the question of the correct treatment of the vertex corrections JTO! . It is 
known that the density response function for the electron gas in the self-consistent GW 
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approximation ^Tj violates the f-sum rule fHlEl- For a given approximation for the self- 
energy the corresponding vertex corrections can be obtained by solving the Bethe-Salpeter 
equation with a specific particle-hole irreducible kernel [201 E|- The actual solution of 
the integral equation for the in-medium vertex is quite complex, a complete solution 
for dressed propagators was given by Kwong and Bonitz |2I] from the solution of non- 
equilibrium Kadanoff-Baym equations in an external field. This procedure is numerically 
quite involved. The initial state is given by a non-equilibrium evolution of the quantum 
transport equations and is not easy to tune. For a quantum well the equations for the 
dressed vertex were solved using equilibrium Green's function technique taking the most 
important part of the particle- hole irreducible kernel [22] • Below we present the first 
solution of the Bethe-Salpeter equation for the dressed vertex in a nuclear system for 
a general momentum and energy of the external field coupled to the density using the 
real-time Green's function formalism. 

Self-consistent approximations for nuclear systems include iterated second Born ap- 
proximation using an effective interaction fl^ El CD] and in medium T- matrix calculations 
using free nucleon-nucleon interactions [211 1221 122] • These schemes use dressed nucleons 
with nontrivial spectral functions and the response function cannot be obtained simply by 
calculating the particle-hole polarization loop. The dressed vertex describing the coupling 
of the external current to the in-medium nucleons is given by a Bethe-Salpeter equation 
(Fig. where K denotes the particle-hole irreducible kernel. For a self-consistent ap- 
proximation to the self-energy one should take for the kernel K the functional derivative 
of the self-energy with respect to the dressed Green's function [201112] K = 6T,/6G. 

In the following we use the self-consistent second Born approximation with the effective 
interaction taken from [22] (Fig. 12]) • The solution of the self-consistent equations for the 
self-energy and dressed propagators in equilibrium is standard 123 HHj and numerical 
easy to implement. In the real time formalism the nucleon retarded (advanced) Green's 
function is expressed by the retarded (advanced) self-energy E^*^") through the Dyson 
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The self-energy is taken in the second direct Born approximation 
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are written using the Fermi distribution f{uj) and the spectral function 
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Equations ((T)), (j2)) and (jD) are solved by iteration with a constraint on the total density 
of the system. In the following we present results for the density p = 0.32fm~^ and the 
temperature of 5MeV. The single-particle width is of about 20 — SOMeV. The polarization 
bubble without vertex corrections is given by 

n<(»(q,fi) = -4z / ^G<>\ci + p,co + n)G>^<\p,co) . (6) 

The retarded and advanced polarization n^*^'^'' can be obtained using a dispersion relation 
from its imaginary part 2Imn'' = 11^ — 11^. 

Let us define the Green's function in a weak external field of momentum q and energy 
Q coupled to the density. For the incoming nucleon momentum p and energy a; the 
outgoing momentum is q + p and the energy uj + fl. We define the smaller (larger) 
Green's functions 

Gf^\q + p,uj + n;p,uj) (7) 
and the retarded (advanced) Green's functions 

G'j^''\q + p,u + n;p,u) , (8) 

with respect to the ordering of the fermion legs in the three-point function Gj, the de- 
pendence on the external energy Q corresponds always to a retarded vertex. They fulfill 
the relation 

G}~G'} = G^ - Gf (9) 

but there is no spectral representation for them. It means that there are three indepen- 
dent Green's functions Gj as expected j2H|- The Green's functions Gj are given by the 
amputated vertex T 

G'i^''\q + p,u + n;p,uj)= G'^''\q + p,uj + fi)r(")(q + p, + fi; p, cu)G"^(")(p, uj) (10) 
and 

Gf''\q + p,u + n;p,u) = G^'(q + P, ^ + f^)r''(q + p, u; + fi; p, u)G<>\p, u) 

+G"(q + p,u + n)T<>\q + p,u + Q;p, cu)G"^(p, u) 
+G<>\q + p,uj + fi)r(q + p, + fi; p, cu)G"(p, cu) . (11) 

The dressed vertex F for the coupling of the external field to the nucleon is the solution of 
the Bethe-Salpeter equation. The particle-hole irreducible kernel for the self-energy given 
by Eq. Q is shown in Fig. |21 

Defining the polarization bubble coupled to the external field 

(^2Tr-^i^ {Gi {q + P + Pi, uj + n + uJi;p + Pi, u + Ui) 
G>«) (pi, cu^) + G<>\p + pi, + cui)G>(<) (pi, cu^; pi -q,uj,- ^]))(12) 



3 



and 



Ifj^''\q + p,uj + n;p,uj) = -Ai j ^^^^^^ (G'^^''^(q + p + pi, + + p + Pi, u; + cJiJ 



(2vr)^ 

G<(»(pi, + G"^("nP + Pi, ^ + ^^^^^(^^(pi, ^i; pi - q, - n) 
+G^^>^(q + p + pi, + + cJi; p + pi, + cJi)G"'(")(pi, cJi) 

+G<(»(p + pi, + cui)G^('^)(pi, pi - q, u;i - fi))(13) 



we can write the Bethe-Salpeter equations for the dressed vertex as 

r<(»(p + q, + fi; p, c.) = ^ / {y\v - Pi)n<(>)(p - pi, ^ - ^0 

^^""(Pi + q, tui + il; Pi, t^i) 
+^(Pi)^(Pi + q)n,<^>^ (pi + q, + fi; pi, ^i)G<(»(p - Pi, ^ - Ur)) (14) 



r^'^) (p + q, a; + 1]; p, ^) = 1 + ^ I (^^^(p - pi)n<(p - pi, a; - 



and 

G'?''''(pi + q, t^i + fi; Pi, cui) 

+^(pi)^(pi + q)n^(pi + q, ^1 + pi, ^^^^"("^(p - pi, - c^i) 

+V\v - - Pi, ^ - ^i)Gi^(Pi + q, ^1 + ^\ Pi, c^i) 

+^(pi)^(pi + q)n?"^(pi + q, + pi, ^i)G'>(p - pi, ^ - ^i)) . (is) 



The equations for the dressed vertex (|T^. (fT3j) and for Green's function in the external 
field ()10|) . (fTTj) are iterated for each given q and f2, using the Green's functions G obtained 
for the correlated system in equilibrium. The irreducible polarization bubble including 
vertex corrections for the density response is calculated from the Green's function Gf 

n[„(q, fi) = -At j ^^G<(p + q, + p, .;) . (16) 

In Fig. m is show the polarization bubble with vertex corrections Tlirr for q = 220MeV. 
The result is very different from the naive one-loop polarization © , the vertex corrections 
are very important. The response function is closer to the response of a noninteracting 
system, but not exactly the same. We have made a calculation of the polarization function 
neglecting the contributions of the terms with 11/ in the Bethe-Salpeter equations (jl4j) . 
(jlSp . This corresponds to neglecting the last two diagrams in the particle- hole kernel 
(Fig. 121). The result is different which shows that the full kernel of the Bethe-Salpeter 
equation must be taken and not only the first term in Fig. El 

Unlike for the function Gj, the real and imaginary parts of the polarization bubble 
U^j.^ fulfill a dispersion relation 

T) ur \ f Imn;,^(q,^i) 

Ren^rr q,^ = - / —— ■ 17 

J TT UJ — UJi + te 
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The dispersion relations between the independently obtained real and imaginary parts of 
the response function constitutes a useful check of the consistency of the approach and 
of the numerics. Another important consistency of the response functions is given by the 
f-sum rule 



The above sum rule is well satisfied for the free response and for the response function 
including vertex corrections, it is severely violated for the naive polarization loop with 
dressed propagators Please note that the sum rule is satisfied also by the solution 
using a simplified kernel in the Bethe-Salpeter equation, although the polarization itself 
is quite different. Our numerical method is quite accurate for q > lOOMeV, where both 
the f-sum rule and the dispersion relation (|17|) are satisfied. Note that in the region where 
the calculations are performed the vertex corrections are still very important, essential 
in guaranteeing the fulfillment of the sum rule. For example at g = 220MeV the naive 
polarization loop with dressed propagators violates the sum rule by a factor 4.9 . 

Adding the Hartree-Fock terms in the self-energy modifies the kernel of the equation 
for the dressed vertex. The Fock term generates an additional interaction line in Fig. |2l 
the exchange term in the random-phase approximation. The Hartree term generates the 
ring series which is summed to give the polarization 



We have checked that adding at the same time a Hartree-Fock term to the self-energy ^ 
and an exchange interaction ladder in the kernel of the Bethe-Salpeter equation conserves 
the sum rule for II^^j,. This sum rule is then conserved after performing the ring summa- 
tion ()19p . although the shape of the response function itself is strongly modified by the 
transformation (fT^ . The irreducible response function with vertex correction is close to 
the noninteracting one, however after the ring summation the results are different due to 
a different tail in Imflj^^ at large energies. The study of the full response with a realistic 
effective Hartree-Fock self-energy will be presented elsewhere. 

We present a solution of the Bethe-Salpeter equation for the vertex of the external field 
coupled to the density in a correlated medium. The approximation chosen for the particle- 
hole interaction is consistent with the nontrivial self-energy used in the Dyson equation. 
It guarantees the fulfillment of the f-sum rule for the response function. To our knowledge 
it is the first solution of the Bethe-Salpeter equation for a general momentum and energy, 
using equilibrium Green's function formalism and the full kernel of the integral equation 
for the dressed vertex. The irreducible polarization with vertex corrections is much closer 
to the Lindhard function than the naive one loop polarization with dressed propagators. It 
is a manifestation of the expected cancellation of the self-energy and vertex corrections. 
The methods here presented can be used to calculate realistic response functions with 
different vertices or to test approximations for the dressed vertices or for the particle-hole 
interaction for the more complicated T-matrix self-energies in nuclear matter. 
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Figure 1: The Bethe-Salpeter equation for the dressed vertex. The particle-hole irre- 
ducible kernel K is denoted by the box and the fat and the small dots denote the dressed 
and the bare vertices respectively. 




Figure 2: The self-energy in the Born approximation. 




Figure 3: The particle-hole irreducible vertex corresponding to the self-energy in Fig. |21 
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Figure 4: The imaginary part of the irreducible polarization as function of the energy 
for q = 220MeV. The solid line denotes the result with vertex corrections, the dotted 
line is the one-loop calculation with dressed propagators (Eq. and the dashed line 
is the polarization for the noninteracting system. The dashed-dotted line denotes the 
polarization with vertex corrections restricted only to the first diagram in Fig. IHl 
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Figure 5: The ratio of the f-sum rule to the expected value as function of the momentum 
in the response function. The lines are the same as in Fig. |3 
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